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Abstract A 

singularly perturbed linear system of second order ordinary differential 
equations of reaction-diffusion type with given boundary conditions is con- 
sidered. The leading term of each equation is multiplied by a small positive 
parameter. These singular perturbation parameters are assumed to be dis- 
tinct. The components of the solution exhibit overlapping layers. Shishkin 
piecewise-uniform meshes are introduced, which are used in conjunction with 
a classical finite difference discretisation, to construct a numerical method 
for solving this problem. It is proved that the numerical approximations ob- 
tained with this method is essentially second order convergent uniformly with 
respect to all of the parameters. 



1 Introduction 

The following two-point boundary value problem is considered for the singu- 
larly perturbed linear system of second order differential equations 

-Eu"(x)+A(x)u(x)=f(x), x € (0, 1), u(0) and u(l) given. (1) 

Here u is a column n — vector, E and A(x) are n X n matrices, 
E = diag(e), e = (ex, ■ ■ ■ , e n ) with < e, < I for all i = 1, . . . , n. The 
£j are assumed to be distinct and, for convenience, to have the ordering 



ex < ■■■ < e n . 
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Cases with some of the parameters coincident are not considered here. 
The problem can also be written in the operator form 

Lu = f , u(0) and u(l) given 

where the operator L is defined by 

d 2 

L = -ED 2 + A(x) and D 2 = —. 

dx z 

For all x e [0, 1] it is assumed that the components aij(x) of A(x) satisfy the 
inequalities 

n 

a-u(x) > /~J a ij( x )\ f° r 1 < i < n, and aij(x) < for i ^ j (2) 

j = l 

and, for some a, 

n 

0<a< min (Va.jfi)). (3) 

a; £[0,1 J ' » 

i<i<„ j=i 

Wherever necessary the required smoothness of the problem data is assumed. 
It is also assumed, without loss of generality, that 

max y/Fi<^— . (4) 

l<z<n 

The norms || V ||= maxi<K„ \Vk\ for any n- vector V, || y \\= sup < 2 ,< 1 \y(x)\ 
for any scalar- valued function y and || y ||= maxi<fe< rl || yu \\ for any vector- 
valued function y are introduced. Throughout the paper C denotes a generic 
positive constant, which is independent of x and of all singular perturbation 
and discretization parameters. Furthermore, inequalities between vectors are 
understood in the componentwise sense. 



For a general introduction to parameter-uniform numerical methods for sin- 
gular perturbation problems, see [I], [2] and [3]. Parameter-uniform numerical 
methods for various special cases of (dJ are examined in, for example, [5] , [5] 
and [7] . For ([T]) itself parameter- uniform numerical methods of first and sec- 
ond order are considered in [8]. However, the present paper differs from [8] 
in two important ways. First of all, the meshes, and hence the numerical 
methods, used are different from those in [8]; the transition points between 
meshes of differing resolution are defined in a similar but different manner. 
The piecewise- uniform Shishkin meshes Mb in the present paper have the ele- 
gant property that they reduce to uniform meshes whenever b = 0. Secondly, 
the proofs given here do not require the use of Green's function techniques, 
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as is the case in [5] . The significance of this is that it is more likely that such 
techniques can be extended in future to problems in higher dimensions and 
to nonlinear problems, than is the case for proofs depending on Green's func- 
tions. It is also satisfying to demonstrate that the methods of proof pioneered 
by G. I. Shishkin can be extended successfully to systems of this kind. 

The plan of the paper is as follows. In the next section both standard and 
novel bounds on the smooth and singular components of the exact solution 
are obtained. The sharp estimates for the singular component in Lemma[7]are 
proved by mathematical induction, while interesting orderings of the points 
established in Lemma [5] In Section 4 piecewise-uniform Shishkin 
meshes are introduced, the discrete problem is defined and the discrete max- 
imum principle and discrete stability properties are established. In Section 
6 an expression for the local truncation error and a standard estimate are 
stated. In Section 7 parameter-uniform estimates for the local truncation er- 
ror of the smooth and singular components are obtained in a sequence of 
theorems. The section culminates with the statement and proof of the essen- 
tially second order parameter-uniform error estimate. 



2 Standard analytical results 

The operator L satisfies the following maximum principle 

Lemma 1. Let A(x) satisfy (0) and f3|). Let ip be any function in the 
domain of L such that ip(0) > and tp(l) > 0. Then ~Lip(x) > for 
all x € (0,1) implies that ip(x) > for all x G [0,1]. 

Proof. Let i*,x* be such that ipi*( x *) — mhij iX ipi(x) and assume that the 
lemma is false. Then ipi* (x*) < . From the hypotheses we have x* $ {0, 1} 
and ?/>■'» 0*) > 0. Thus 

n 

i=i 

which contradicts the assumption and proves the result for L. ■ 

Let A(x) be any principal sub-matrix of A{x) and L the corresponding 
operator. To see that any L satisfies the same maximum principle as L, 
it suffices to observe that the elements of A(x) satisfy a fortiori the same 
inequalities as those of A(x). 

Lemma 2. Let A(x) satisfy (0) and $3$. If xp is any function in the domain 
of L, then for each i, 1 < i < n, 

\i>i{x)\ < max (|| V(0) ||,|| </>(!) II, - II W \\) , a €[0,1]. 
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Proof. Define the two functions 

9 ± (x) = maxjll tf(0) ||, || V(l) II, \ II W ||}e ± ^(x) 

where e = (1, 1) T is the unit column vector. Using the properties of 
A it is not hard to verify that #±(0) > 0, 6 ± (l) > and L6 ± (x) > 0. 
It follows from Lemma Q] that G ± (x) > for all x £ [0, 1]. ■ 

A standard estimate of the exact solution and its derivatives is contained in 
the following lemma. 

Lemma 3. Let A{x) satisfy (0) and (Q) and let u be the exact solution of 
(QP. Then, for each i = 1 . . . n, all x £ [0, 1] and k = 0, 1, 2, 



K (fc) (x)|<C £ ^(||u(0)|| + ||u(l)|| + ||f||) 
\uf\x)\ < Cfer*(|| u (0)|| + ||u(l)|| + ||f || + ^l\\f 



\uf\x)\ < Cer2(|| u ( )|| + || U (1)|| + ||f|| + £i ||f"||). 

Proof. The bound on u is an immediate consequence of Lemma [5] and the 
differential equation. 

To bound u'^x), for all i and any x, consider an interval N x = [a, a + y^i] 
such that x £ N x . Then, by the mean value theorem, for some y £ N x , 

,, , Ui(a + J~el) - Ui(a) 

u i(y) = —}= 

and it follows that 

|<( 2 /)|<2 £ p|| Ui ||. 

Now 

pX pX 

u'(x) = u'(y) + / u"(s)ds = u'(y) + E' 1 / (-f(s) + A(s)u(s))ds 
Jy Jy 

and so 

K(*)| < \<{y)\ + CezXWM + ||u||) f ds < Ce;H\\m + ||u||) 



from which the required bound follows. 

Rewriting and differentiating the differential equation gives u" = E~ 1 (Au — 
f), u< 3 ) =E- 1 (Au' + A'u-f), uW = E-^Au" +2A'u' + A"u-f"), and 
the bounds on u'/, u[ 3 \ u\ follow. ■ 

The reduced solution Uo of ((TJ) is the solution of the reduced equation Avlq = 
f. The Shishkin decomposition of the exact solution u of JT]) is u = 
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v + w where the smooth component v is the solution of Lv = f with 
v(0) = Uq(O) and v(l) = u (l) and the singular component w is the 
solution of Lw = with w(0) = u(0) - v(0) and w(l) = u(l) - v(l). 
For convenience the left and right boundary layers of w are separated using 
the further decomposition w = w' + w r where Lw' = 0, w'(0) = u(0) - 
v(0), w'(l) = and Lw r = 0, w r (0) = 0, w r (l) = u(l) - v(l). 
Bounds on the smooth component and its derivatives are contained in 

Lemma 4. Let A(x) satisfy p|) and (0). Then the smooth component v and 
its derivatives satisfy, for all x £ [0,1], i = 1, . . . n and k = 0, ... 4, 

\v i t\x)\<C{l+e\- ,i ). 

Proof. The bound on v is an immediate consequence of the defining equations 
for v and Lemma [2] 

The bounds on v' and v" are found as follows. Differentiating twice the 
equation for v, it is not hard to see that v" satisfies 

Lv" = g, where g = f " - A"\ - 2A'v h (5) 

Also the defining equations for v yield at x = 0, x = 1 

v"(0) = 0, v"(l) = 0. (6) 

Applying Lemma [2] to v" then gives 

||v"||<C(l + ||v'||). (7) 

Choosing i* , x* , such that 1 < i* < n, x* € (0, 1) and 

v' i ,(x*) = \\v'\\ (8) 

and using a Taylor expansion it follows that, for any y € [0, 1 — x*] and some 
T], x* < i] < x* + y, 

(x* +y)= v t . (x*)+y v[, (x*) + V - v'>, fa). (9) 
Rearranging © yields 

) = " 9 U i*W ( 10 ) 

y L 



and so, from flU) and ([10 



|v'||<^||v|| + f||v"||. (11) 

y L 



Using (|TT|) . (J7]) and the bound on v yields 
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(l-^)||v"||<C(l + -). (12) 

2 y 

Choosing y = min(^,l — x*), (IT2l) then gives ||v"|| < C and (fTTj) gives 
||v'|| < C as required. The bounds on v( 3 ),v' 4 ) are obtained by a similar 
argument. ■ 



3 Improved estimates 

The layer functions B\ 1 B\, Bi, i = 1, . . . , n, , associated with the solution 
u, arc defined on [0, 1] by 

B\(x) = e- x ^\ B\{x) = B\{1 - x), B t (x) = B\{x) + BV(x). 

The following elementary properties of these layer functions, for all 1 < i < 
j < n and < x < y < 1, should be noted: 

(a) Bj(a;) < B\{x), B\{x) > B\{y), < B\{x) < 1. 

(b) BT(x) < Br(x), BV(x) < BT(y), < Bl(x) < 1. 

(c) Bi(x) is monotone decreasing (increasing) for increasing x £ [0, ^]([^, 1])- 

(d) Bi{x) < 2B\(x) for x G [0, §]. 

Definition 1. For £>,-, B^, each i,j, 1 < i ^ j < n and each s, s > 0, the 

B i M:h b\(x { :}) 



point x[ S j is defined by 



It is remarked that 



5[(l-xg) SJ(l-xg) 



£ ■ 



(13) 



(14) 



In the next lemma the existence and uniqueness of the points x\ s ^ are shown. 
Various properties are also established. 

Lemma 5. For all i,j, such that 1 < i < j < n and < s < 3/2, the points 
Xij exist, are uniquely defined and satisfy the following inequalities 

5M>5jM , e|Mg) , £M < 5M, l6(lS , 11 . (15) 

i 3 i 3 

Moreover 

41 < if * + K j and zg < zg +1 , if i < j. (16) 
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Also 

zg <2s^ and zg e (0, i) if i < j. CI 7 ) 

Analogous results hold for the B\, Z?J and t/ie points 1 — ^g- 



Proof. Existence, uniqueness and (| 15[) follow from the observation that the 
ratio of the two sides of (TT51) . namely 



£ j £ j , r- , 1 1 ^ 

iTBiw = < exp( - vSl< vS-^J )) ' 

is monotonically decreasing from the value |£ > 1 as x increases from 0. 

The point a;g is the unique point x at which this ratio has the value ', 
Rearranging (|13[) , and using the inequality In x < x — 1 for all x > 1 , gives 



*8 = 2s 



2s ln(^) /i~ 
-^—<2aJ^-, (18) 



which is the first part of (fTTj) . The second part follows immediately from this 
and flU). 

To prove (|16[) . writing ^/e^ = exp(— pk), for some p^ > and all k, it follows 
that 

^(s) = 2s{pi-pj) 
hj ^/a(exppi - exppj) 

(s) Is) 

The inequality x\ ! < ■ is equivalent to 

Pt - Pj Pi+i- Pj 



exppi - exppj exppi + i - expp,, 

which can be written in the form 

(Pi+i - P 3 ) exp(pi -pj) + (pi - Pi+i) - (pi - Pj) exp(p J+ i - Pj) > 0. 

With a — pi — pj and b — pi+i — pj it is not hard to see that a > b > and 
a — b = pi — Pi+i. Moreover, the previous inequality is then equivalent to 

exp a — 1 exp 6 — 1 
a b 

which is true because a > b and proves the first part of (| 16[) . The second part 
is proved by a similar argument. 
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The analogous results for the B r - and the points 1 — x\J are proved by 
a similar argument. ■ 

In the following lemma sharper estimates of the smooth component are pre- 
sented. 

Lemma 6. Let A(x) satisfy JJ|) and |3|). Then the smooth component v 
of the solution u of ([!]) satisfies for i = 1, • • • , n, k = 0, 1, 2, 3 and x G fi 



\v\ h \x)\ < C 1 + E 



B q {x) 



q=l Cq 

Proof. Define a barrier function 

^(x) = C[l + B n (x)]e ± v (fc) (x), fc = 0,l,2 and x G 77. 

Using Lemma [D we find that LV ± (a;) > and t/> ± (0) > 0, V ± (l) > 

for proper choices of the constant C. 

Thus using Lemma 0] we conclude that for k = 0, 1, 2, 



I«i*°(a0l < C[l + S„(a:)], i G /?. 



(19) 



Consider the system of equations ([5]), ([5]) satisfied by v", and note that 

II g' II — C from LemmaHl 

For convenience let p denote v" then 

LP = g, P(0) = 0, p(l) = 0. (20) 

Let q and r be the smooth and singular components of p given by 

Lq = g, q(0) = A(0)- X g(0), q(l) - Afl)- 1 ^) 

and 

Lr = 0, r(0) = -q(0), r(l) = -q(l). 
Using Lemmas S] and [7] we have, for i = 1, • •• ,n and x G £2, 

Wi{x)\ < c, 
\r'i(x)\ < C 



Bi(x) 



B n (x) 



Hence, for x G fl and « = !,••■ , n, 



\v'/'(x)\ < \p' t (x)\<C 



Bi{x) 



B n {x) 



(21) 



From dHj) and (|2Tj), we find that for k = 0, 1, 2, 3 and x G (2, 
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\vi k) (x)\ < C 1 + 4 *Bi(x) + 



T B n {x) 



Remark : It is interesting to note that the above estimate reduces to the 
estimate of the smooth component of the solution of the scalar problem given 
in [I] when n = 1. 

Bounds on the singular components w ; , w r of u and their derivatives are 
contained in 

Lemma 7. Let A(x) satisfy (0) and (dp. Then there exists a constant C, such 
that, for each x G [0, 1] and i = 1, . . . , n, 



\w\{x)\ < CB l n (x), w l /{x) < CY, 



n nl 



B l q (x) 



1,11 I N 

W:' (x) 



: C 



B\{x) 



q=i 



sat' 



(,:) 



1,(3) 



(x) 



< 



3/2 



< 



9=1 



B\{x) 



Analogous results hold for and its derivatives. 

Proof. First we obtain the bound on w l . We define the two functions 9 ± = 
CB l n e ± w*. Then clearly 9 ± (0) > 0, 9 ± (1) > and L9 ± = CL(B l n e). 
Then, for i = 1, . . . , n, (L0 ± ) 4 = C(£™ =1 a 4ji - af^)B l n > 0. By Lemma [TJ 
± > 0, which leads to the required bound on w'. 

Assuming, for the moment, the bounds on the first and second derivatives 
w/ and wS' , the system of differential equations satisfied by w' is differen- 
tiated twice to get 



-Ew l 'W +Aw l >" + 2A\ 



0. 



The required bounds on the w 1 ^^ follow from those on w[, w 1 / and w\'" . It 
remains therefore to establish the bounds on w 1 / ,w\'" and w/" , for which 
the following mathematical induction argument is used. It is assumed that 
the bounds hold for all systems up to order n — 1. It is then shown that the 
bounds hold for order n. The induction argument is completed by observing 
that the bounds for the scalar case n = 1 are proved in pQ. 

It is now shown that under the induction hypothesis the required bounds 
hold for w 1 / ,w\'" and w/" . The bounds when i = n are established first. The 
differential equation for w l n gives e n w 1 ^' = (4w )„ and the required bound 



on w 1 '" follows at once from that for w 1 . For 



it is seen from the bounds 

l,i> 



in Lemma El applied to the system satisfied by w ( , that [x)\ < Ce i 2 . 
In particular, |w^'(0)| < Ce n 2 and |u^'(l)| < Ce n 2 . It is also not hard 
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to verify that Lw ! '' = — A'w. Using these results, the inequalities £j < 
e n , i < n, and the properties of A, it follows that the two barrier functions 
± = CE-iB l n e±w l >' satisfy the inequalities 9 ± (0) > 0, ± (l) > and 

L0 ± > 0. It follows from Lemma Q] that 9 ± > and in particular that its 

_ i 

n th component satisfies < Ce„ 2 B l n (x) as required. 

Now, consider 

-EnW l r i / '(x)+a n i(x)w{(x)+a n 2(x)w l 2 (x)-\ \-a nn (x)w l n (x) = f n (x). (22) 

Differentiating (f22|) once, we get 



l wll ( - 3 \x) = f' n {x) - V (a nj (x)w l j (x)) 



|^ (3) (x)| 



< Ce- 1 



< Ce: 1 



B[(x) 



B l n {x) 



3/2 

q=l £ q 



To bound u;.'', w/' 



and w' ,< ' 3 ' ) for 1 < i < n — 1 introduce w 1 = 
(u^, . . . , tf*,^). Then, taking the first n — 1 equations satisfied by w', it 
follows that 

Aw' = g, 



-Ew 1 '" 



where E, A is the matrix obtained by deleting the last row and column 
from E, A, respectively, and the components of g are gi = —ai, n w l n for 



1 < i < n— 1. Using the bounds already obtained for w n , w. 



i.i 



and w 1 '"' 



it is seen that g is bounded by CB l n (x), g' by C 



g" by C 



and g" 



by C 



Q ^ *^ ' 11 

— ^2~. The boundary conditions for w are w (0) = u(0) — u (0), 

9=1 e g 

w z (l) = 0, where u° is the solution of the reduced problem u° = A~ 1 f 1 and 
are bounded by C(|| u(0) || + || f(0) ||) and C(|| u(l) || + || f(l) ||). Now 
decompose w' into smooth and singular components to get 



w' = q 



w''' = q' 



Applying Lemma [T] to q and using the bounds on the inhomogeneous term g 
and its derivatives g', g" and g( 3 ) it follows that |q'(x)| < C Bn S^ , |q"(a^)| < 



Numerical solution of a reaction-diffusion system 11 

(jRjAsl anc i |q'"(a;)| < — \pT' Using mathematical induction, assume 

9=1 £ 9 

that the result holds for all systems with n — 1 equations. Then Lemma [7] 
applies to r and so, for i = 1, . . . , n — 1, 

M(*)l<c£^p. K'(,)|<c^^, |rf(*)|<c£^J 



9 =i V "« ?= i q=l £9 

Combining the bounds for the derivatives of qi and r^, it follows that 

k^)i<cE^, \ w y{x)\<c±^ \ w r'( X )\<c± B 'i {x 



F ■ . . - - ■ - 3/2 • 

A/ -4 - c 9 1 F« 

9=2 V * q—i 1 q—\ ^9 

Thus, the bounds on w 1 / , w/' and w 1 /" hold for a system with n equations, 
as required. A similar proof of the analogous results for the right boundary 
layer functions holds. ■ 



4 The Shishkin mesh 



A piecewise uniform mesh with N mesh-intervals and mesh-points {xi}fL 
is now constructed by dividing the interval [0, 1] into 2n + 1 sub- intervals as 
follows 

[0, n] U • • • U (r»_i,r n ] U (r„, 1 - r„] U (1 - r„, 1 - r„_ x ] U • • • U (1 - n, 1]. 

The n parameters r^, which determine the points separating the uniform 
meshes, are defined by 

1 2.,/^lniv} (23) 



and for k = 1, 



4' 



T fc = mini^±i,2 A /-ln7V!>. (24) 
Clearly 

1 3 

< n < .. . < r„ < -, - < 1 - r„ < . . . < 1 - Tl < 1. 

Then, on the sub-interval (r„,l— r„] a uniform mesh with y mesh-intervals 
is placed, on each of the sub-intervals (rfe, Tfc+i] and (1 — t^+i, 1 — Tfe], fc = 
1, . . . , n — 1, a uniform mesh of „„^ +2 mesh- intervals is placed and on 
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both of the sub-intervals [0,ti] and (1 — ri,l] a uniform mesh of t^tt 
mesh-intervals is placed. In practice it is convenient to take 

N = 2 n+p+1 (25) 

for some natural number p. It follows that in the sub- interval [r^-i, Tk] there 
are N/2 n ~ k+3 = 2 k+p ~ 2 mesh-intervals. This construction leads to a class of 
2™ piecewise uniform Shishkin meshes Mb, where b denotes an n-vector with 
bj = if Ti = and 6^ = 1 otherwise. From the above construction it clear 
that the only points at which the meshsize can change are in a subset Jb of 
the set of transition points Tb = {Tfe}^ =1 U {1 - t^Lj. It is not hard to see 
that the change in the meshsize at each point Tk is 2 n ~ k+3 (dk — dk-i), where 
dk = -^2^ ~ T k f° r 1 < k < n, with the conventions do — 0, T n +i = 1/2. 
Notice that dk > and that bk = if and only if dk — 0. It follows that Afb 
is a classical uniform mesh when b = 0. 
The following notation is now introduced: Hj 

Xj-i, Sj = xj + i — Xj-i, Jb = {xj : Hj — hj ^ 0}. Clearly, Jb is the set 
of points at which the meshsize changes and Jb C Tb. Note that, in gen- 
eral, Jb is a proper subset of Tb. Moreover, if bk = then Hk < hk and if 
bk = bk-i = then Hk — hk- In the latter case, it follows that the meshsize 
does not change at Tk or 1 — Tk. 
It is not hard to see also that 

T fc <C\/efchi7V, 1 < k < n, (26) 

h k = 2 n - k+3 N~\T k - t*_i), H k = 2 n - k+2 N- 1 ( Tk+1 - T k ), (27) 
S j =H j + h j <Cmax.{H j ,h j }, 1 < j < N - 1, (28) 
T k = 2- (] - k+1) T j+l when b k = ■ ■ ■ = bj ■ = 0, 1 < k < j < n (29) 

and 

B l k(T k ) = B r k (l - T k ) - N- 2 when b k - 1. (30) 
The geometrical results in the following lemma are used later. 
Lemma 8. Assume that bk = 1. Then the following inequalities hold 

Xk-i.k - T k - hk for 1 < k < n. (31) 

B\{r k ) < J_ for x < , < n (32) 



B\{T k - h k ) < CB l q (Tk) for 1 < k < q < n. (33) 

Proof. To verify (|3T|) note that by Lemma [5] 

(s) 9 Vgfe _ gTfc _ sT k n 

Xk - 1 ' k _ IniV ~ (n + p + l)ln2 " 2 ' 
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Also, 

AT ~ 2 2 

( s) 

It follows that a^-i k + h k < r k as required. 

To verify (j3"2")l note that if i > k the result is trivial. On the other hand, if 
i < k, by (f3Tj) and Lemma El 

B\{r k ) < ^ gUjg) < 1 



Finally, to verify (f33|) note that 

^ = (n -r fc _i)2"- fe+3 iV- 1 < r fc 2 n - fe+3 iV- 1 = A /^2 rl - fe+4 iV- 1 lnAr. 

V a 



and 



e 2»-*+'jv-M„Jv = (JV * )2 "-*+« < Cj 



so 

A^ftfe < /1^2"- fc+4 7V^ 1 ln7V < 2 n - fe+4 iV- 1 lniV < C 

V ^9 V ^9 

since k < q. It follows that 

Bjfa - &k) - B l q (T k )eV^ hk < CB l q {r k ). 

as required. ■ 



5 The discrete problem 

In this section a classical finite difference operator with an appropriate 
Shishkin mesh is used to construct a numerical method for (|T|), which is 
shown later to be essentially second order parameter-uniform. In the scalar 
case, when n = 1, this result is well known. In [7] it is established for general 
values of n in the special case where all of the singular perturbation parame- 
ters are equal. For the general case considered here, the error analysis is based 
on an extension of the techniques employed in [3J. It is assumed henceforth 
that the problem data satisfy whatever smoothness conditions are required. 
The discrete two-point boundary value problem is now defined on any mesh 
Mb by the finite difference method 



-ES 2 V + A{x)U = f(x), U(0)=u(0), U(l) = u(l). (34) 
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This is used to compute numerical approximations to the exact solution of 
(HJ. Note that ([M]) can also be written in the operator form 

L N U = f, U(0) = u(0), U(l) = u(l) 

where 

L N = -E5 2 + A(x) 
and 5 2 , D + and D~ are the difference operators 

D+U(xj) - D-\J(xj) 



, , N U(x,-+i) -U(xi) , s Ufa;,-) - U(x,-_i) 

D + U(xj) = V 3+1 ' ^ and D-\J{ Xj ) = 3 ' J ■ 

with /ij = J — J+1 , hj = Xj — Xj-\. 

The following discrete results are analogous to those for the continuous case. 

Lemma 9. Let A(x) satisfy (0) and Then, for any mesh function the 
inequalities *(0) > 0, > and L N *( Xj ) > for 1 < j < N-l, 

imply that > for < j < N. 

Proof. Let i*,j* be such that \Pi*(xj*) = minjj &i(xj) and assume that the 
lemma is false. Then ^* (xj* ) < . From the hypotheses we have j* ^ 0, N 
and &i.(xj.)- &i*(xj.-i) < 0, &i*(xj* +1 ) - &i»{x r ) > 0, so 5 2 &u(xj*) > 
0. It follows that 

n 

(L N i&(xj*)) i:t = -£i*5 2 &„( X j*) + ^2 ai *> k(Xj*)&k(%j*) < 0, 

fe=l 

which is a contradiction, as required. ■ 

An immediate consequence of this is the following discrete stability result. 

Lemma 10. Let A(x) satisfy and (0). Then, for any mesh function \&, 

|| 9( Xj ) || < max{||¥(0)|W|¥(l)||,~||L*¥||}, < j < N. 

Proof. Define the two functions 

®±(x d ) = max{||¥(0)||, ||*(1)||, -||L N *||}e± *( Xj ) 

a 

where e = (1, ... , 1) is the unit vector. Using the properties of A it is not 
hard to verify that ©±(0) > 0, ± (1) > and L N ± (x J ) > 0. It follows 
from Lemma H that ® ± (x 3 ) > for all < j < N. . 
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The following comparison result will be used in the proof of the error estimate. 



Lemma 11. Assume that the mesh functions 4? and Z satisfy, for j = 
1 ... N — 1, 

||Z(0)||<#(0), ||Z(1)|| <*(!), \\(L N )(Z( Xj ))\\<(L N )mxj))- 
Then, for j = ... N , 

||Z(x,)||e<*(x,). 
Proof. Define the two mesh functions \J ,± by 

* ± = *±Z. 

Then * ± satisfies, for j = 1 ... N - 1, 

* ± (0) = * ± (1)=0, (L) A '(* ± )(x J ) > 0. 
The result follows from an application of Lemma [9] ■ 

6 The local truncation error 

From Lemma [POl it is seen that in order to bound the error ||U — u|| it suffices 
to bound h N (XJ — u). But this expression satisfies 

L N (XJ - u) = L W (U) - L w (u) = f L N (u) = L(u) - L w (u) 

= (L - L N )u = -E(5 2 - D 2 )u 

which is the local truncation of the second derivative. Let V, W be the dis- 
crete analogues of v,w respectively. Then, similarly, 

L N (V - v) = -E(5 2 - £> 2 )v, L N (W - w) = -E{8 2 - £> 2 )w. 

By the triangle inequality, 

|| L W (U - u) || < || L N (V v) || + || L N (W w) || . (35) 

Thus, the smooth and singular components of the local truncation error can 
be treated separately. In view of this it is noted that, for any smooth function 
ip, the following three distinct estimates of the local truncation error of its 
second derivative hold: 
for Xj €E Mb 

\(S 2 - D 2 )ip(xj)\ < Cmax\ip"(s)\, (36) 

s 6 Ij 

and 
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\(S 2 - D 2 )ij(x 3 )\ < CJ,-max|V (3) (*)|. (37) 

for Xj Jb 

\(5 2 - D 2 U(x,)\ < CS 2 

sei 

for Tfe E J b 



\(6 2 -D 2 )^(x,)\ < C^max|V (4) (s)|, (38) 



\(5 2 - D 2 ^(T k )\ < C{\H k -h k \.\^\ Tk )\+5 2 m^\^\ S )\). (39) 



7 Error estimate 

The proof of the error estimate is broken into two parts. In the first a theorem 
concerning the smooth part of the error is proved. Then the singular part of 
the error is considered. A barrier function is now constructed, which is used 
in both parts of the proof. 

For each k € lb, introduce the piecewise linear polynomial 

x 

— , < X < Tfe. 
Tfe 

e k ( X ) ={ 1, Tfe < X < 1 - Tfe. 
1 — X 

, l-T k <X<l. 

Tfe 



It is not hard to verify that, for each k e lb, 



2e, 

L N {e k ( Xj )e)i > { ~ ' T k (Hk~+h k ) 



N(n ( \ \ ^ ) a ^ 777 — nTT' if ^ = T k ^ Jb 



a9 k (xj), if ^ J b . 

On the Shishkin mesh Mb define the barrier function $ by 

*(xj) = CN- 2 (\nN) 3 [l + Ok(xj)}e, (40) 

fce/ b 

where C is any sufficiently large constant. 
Then $ satisfies 

< $i(xj) < CN~ 2 (\nN) 3 , l<i<n. (41) 
Also, , for Xj £ Jb, 

{L N <f>( Xj ))i > CN- 2 (\nN) 3 (42) 

and, for r k e Jb, 
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(L N *&(Tk))i > C(l + —J±——)(N- 1 In iV) 2 , 



fe^{H h + hk ) 
from which it follows that, for r k 6 Jb and -fffc > /ifc, 

(L N ^(r k ))i >C(N- 2 + — ^ = N~ 1 ]nN) (43) 



and, for r k e Jb and H k < h k , 

{L N ${T k )) l >C{N- 2 + — N- x \nN). (44) 

£k 

The following theorem gives the error estimate for the smooth component. 

Theorem 1. Let A{x) satisfy and (0). Let v denote the smooth com- 
ponent of the exact solution from (QP and V the smooth component of the 
discrete solution from \3J$ . Then 

||V - v|| < CiV- 2 (lniV) 3 . (45) 

Proof. An application of Lemma[TT]is made, using the above barrier function. 
To prove the theorem it suffices to show that the ratio 



R(vi(xj)) = — TTFW^j~Trr ' ' x j e M *> 

satishcs 



R(vi(xj)) < C. (46) 

For Xj Jb the bound (|46| follows immediately from Lemma IH (|38|) and 
(ESD- 

Now assume that Xj = e 4- The required estimates of the denominator of 
R{vi{Tk)) are (|43|) and (l44|) . The numerator is bounded above using Lemma 
[6]and (|37[) . The cases b k = 1 and b k = are treated separately and the 
inequalities ([26]l. (|27| . (p5]> . ([30]) and <(33]) are used systematically. 
Suppose first that bk = 1, then there are four possible subcases: 

i < k, H k > hk, R(vi(T k )) < Ce k+1 . 

H k < h k , R{v t {T k )) < Ce k . 
i>k, H k > h k , R(vi(r k )) < Ce k+H f§. (47) 

H k < h k , R(vi(T k )) < Ce k 

Secondly, if 6fc = 0, then b k -\ = 1, because otherwise r k ^ Jb, and further- 
more H k < h k . There are two possible subcases: 

i < k - 1, H k < h k , RMn)) < Ce k {^= + 1). 
* > k - 1, H k <h k , R(v t (n)) < Ce k , ^ (4M 
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In all six subcases, because of the ordering of the £j, it is clear that condition 
P5|) is fulfilled. This concludes the proof. ■ 

Before the singular part of the error is estimated the following lemmas are 
established. 

Lemma 12. Let A(x) satisfy $Bj) and |[3J). Then, on each mesh Mb, for 1 < 
i < n and 1 < j < N , the following estimates hold 

S 2 

\si(S 2 — D 2 )w\(x j)\ < C^- for Xj£J h . (49) 

£i 

An analogous result holds for the . 

Proof. When Xj ^ Jb, from ((38]) and Lemma [3 it follows that 



\ £z{ 5 2 - D 2 )w[{ X] )\ < CS 2 max \e z w l ^\s)\ 

J s e i j 

" Bi{s) CS 2 
< CS 2 max V < — °- 



as required. ■ 

In what follows fourth degree polynomials of the form 

\^ {x-x e ) k u{k) 

Pi-fi{x) = 2^ — ^ — w i w 

fe=0 

are used, where 9 denotes a pair of integers separated by a comma. 

Lemma 13. Let A(x) satisfy (0) and |[3J) and assume that Mb is such that 
b k = 1 for some k, 1 < k < n — 1. Then, for each i,j,l<i<n, l<j<N 
there exists a decomposition 

k+l 
9=1 



for which the following estimates hold for each q and r,l<q<k,0<r<2, 
\e iW ^ 2 \ Xj )\<Ce-iB l q { Xj ) 



and 



\s.w% +1 (x 3 )\<C £ MJ +1 (x,)|<C x: 

q=k+l V « 9 =fc+l 9 

Furthermore, for X j (£ Jb, 
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\s % (S 2 ~ D 2 ) W \( Xj )\ < C(B l k ( Xj ^) + A) (50) 

and, for r k e J b , 

M<5 2 - D 2 )w\{r k )\ < C{ B l k (r k h k ) + -J^). (51) 

Analogous results hold for the w\ and their derivatives. 
Proof. Consider the decomposition 

fe+l 



E 

m— 1 



where the components are defined by 



Wik+1 = J Pi;k,k+l on [0,4^+1 ) 
I w\ otherwise 



and for each m, k > m > 2, 

10. 



Pi;m-l,m On [0,X^i lm ) 

' - Egim+i otherwise 



and 

fc+i 

= - KJ; )9 on [0, 1]. 

q=2 

From the above definitions it follows that, for each m, 1 < m < k, Wi 

To establish the bounds on the fourth derivatives it is seen that: 
for x € [^fcfc+i' 1]) Lemma [7] and x > x k X \ +1 imply that 



MliMI = \e lW ^\x)\ < C ]T ^ < C £ 

9=1 £q q=k+l ° q 



1 B\[x) ^ „ ^ B' q (x) 



for x G [0,x k k+1 ]. Lemma [7] and x < x k ' k+1 imply that 

9=1 £ « 9 =fc+l £ « 9 =fc+l £ « 

and for each m = k, ... , 2, it follows that 
for x € K i^fm = 0; 
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for x € [^m-i.m'^m.m+l]' Lemma [7J implies that 

i^wi<i^ (4) wi+ E mSwi< c 'E £ 7^^ c '% m ; 

q—m+1 q=i 

for x € [0, a;^,— l m]> Lemma[7jand a; < aij^Li m imply that 

l (4)/ \ I I 1,(4,), (1) ^, , ^(^rn-l.m) ^ ^,^m( x rn-l,m) ^ r ,B 

q=l E 1 £,n 

for xe [45,1], =0; 

for x € [0, a; 12]) Lemma [7] implies that 

ie,«,g(*)i < i^' w (*)i+Ei^£!(«)i ^ C E^ < 

For the bounds on the second and third derivatives note that, for each m, 
1 < m < k : 

for x € [x ( ^] m+1 , 1], = = ujjfj,; 

for x G [0, a^ m+1 ], / J(s)<fe = e 4 wg>^ m+ J-e^™^) = 

J X 

(3) / \ 

and so 



|e<«fi,(aOI< / ' | £iW W( S )|d S <— / ' S^(*)ds<C 

cm Jx 



In a similar way, it can be shown that 

\e iW 'l m (x)\<CB l m (x). 
Using the above decomposition yields 

fc 

\ £l (S 2 - D 2 )w\{x 3 )\ < E h(S 2 - D 2 )wi <q (xj)\ + \ £l {5 2 - D 2 )w itk+1 ( Xj )\. 

q=l 

For Xj $l Jb, applying ((55)) to the last term and to all other terms on 
the right hand side, it follows that 

k 

\e t (S 2 -D 2 )w\(x 3 )\ < C(Emax| £ i< g (s)l + ^ max le^J+i («)!)• 
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Then (|50[) is obtained by using the bounds on the derivatives obtained in the 
first part of the lemma. 

On the other hand, for Xj = tj. G J^, applying ([57)1 to the last term and (|36p 
to the other terms, (|51[) is obtained by a similar argument. The proof for the 
w\ and their derivatives is similar. ■ 

In what follows third degree polynomials of the form 



are used, where 9 denotes a pair of integers separated by a comma. 

Lemma 14. Let A(x) satisfy (0) and |[3J) and assume that Mb is such that 
bk = 1 for some k, 1 < k < n— 1. Then, for each i, j, 1 < i < n, 1 < j < N 
there exists a decomposition 




k+l 




m—l 



for which the following estimates hold for each m, 1 < m < k, 



Bi 




B> 




e. 



rn 



and 




Furthermore 




(52) 



Analogous results hold for the w\ and their derivatives. 

Proof. The proof is similar to that of Lemma fT3l with the points x^J replaced 
by the points x\ ! . Consider the decomposition 



fe+i 




m— 1 



where the components are defined by 
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and for each m, k > m > 2, 

l',n,,-l,m 011 [Mm-l,m) 

_ fc+1 

w\ — Wi iq otherwise 

q—m+l 

and 

fc+i 

Wi,i = w\ - ^ Wj,q on [0,1]. 

g=2 

From the above definitions it follows that, for each m, 1 < m < k, Wi 
on [xg^.l]. 

To establish the bounds on the third derivatives it is seen that: 
for x € [^fefc+i) 1]) Lemma [7] and x > %kk+i i m Ply that 



\ w % +1 ( X )\=\^ 3 \ X )\<cj:^l<c y, 



B l q {x) ^ n " B' q (x)_ 



3/2 — Zl^ 3/2 

q=k+l £ Q 



for x G [0> x iffe+i]> Lemma [7] and x < x^^ t imply that 

i«SU«)i = ki^oOi < E ■ E L1 ~r^- E 



n R i ^ (3/2) n n nl (A^/ 2 ) s „ R ; , . 

i, (3)^(3/2) m ^ -"gV^fc.fc+li J g^fc,fc+li a q\ x ) 

3/2 - 2-^, 3/2 - 2-~i 3/2 ' 

8=1 £ 1 q=k+l £ Q q=k+l £ 1 



and for each m = k, ... , 2, it follows that 

forxe [a^+i.l], ^S = 0; 

for x € [x^ 2 ^ m , Xm^m+i\i Lemma [Jj implies that 



for x € [0?^m- 2 im]> Lemma [3 and x < E^-im imply that 

= k (3) («,ji < ^e q v < g m( 3 m /2 M ^ 



for x € [x^ 2) ,l], to^ = 0; 

for x £ [0, ]i Lemma [Jj implies that 
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^+1 n TdI ( \ T)l ( \ 

\w${x)\ < \w^( X )\ + 5>£>(x)| <C^^< 

q=2 q—1 £ q £ 1 

For the bounds on the second derivatives note that, for each m, 1 < m < k 
forxe w'l m = 0; 

for x G [0,x^ +1 ], f^w<%(s)da = < m (*^ + i) " <*» = 

-<m( a; ) 

and so 

K m (x)\ < f m ' m+1 \w^(s)\ds < -gj f "'"^ i?L( S )d S < C%^. 
Finally, since 

k 

\ £l (S 2 - D 2 )w\{ Xj )\ < |£i(5 2 - D 2 )w hm {x,)\ + \ £l {6 2 - D 2 )w itk+1 (xj)\, 



using (|37|) on the last term and (|36|) on all other terms on the right hand 
side, it follows that 



\e t (S 2 - D 2 )w\{ X] )\ <C{J2 max h< m ( S )| + ^ max | £< i«g +1 (*)!)• 

m— 1 

The desired result follows by applying the bounds on the derivatives obtained 
in the first part of the lemma. The proof for the w\ and their derivatives is 
similar. ■ 

Lemma 15. Let A( X ) satisfy £|P and |[3J). Then, on each mesh Alb, the fol- 
lowing estimate holds for i = 1, . . . , n and each j ' = 1, . . . , iV, 

\£i(8 2 — D 2 )wl( X j)\ < CB l n { X j_x). 

An analogous result holds for the w\. 

Proof. From ([5B]) and Lemma for each i = 1 , . . . , n and j = 1 , . . . , N, 
it follows that 

\ei{8 2 - D 2 )w\( Xj )\ < C max\siwl'"(s)\ 

q—'i ^ 

The proof for the w\ and their derivatives is similar. ■ 
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The following theorem provides the error estimate for the singular compo- 
nent. 

Theorem 2. Let A(x) satisfy fl[) and Let w denote the singular com- 
ponent of the exact solution from (QJ) and W the singular component of the 
discrete solution from J3^| ). Then 



Proof. Since w = w + w r , it suffices to prove the result for w and w r sep- 
arately. Here it is proved for w' by an application of Lemma 111! A similar 
proof holds for w r . 
The proof is in two parts. 

First assume that Xj ^ Jb. Each open subinterval (Tfc,Tfc + i) is treated sepa- 
rately. 

First, consider Xj £ (0,ti). Then, on each mesh Mb, 6j < CN^ 1 ^ and the 
result follows from (|26p and Lemma ll"2"1 

Secondly, consider Xj € (ri,r 2 ), then n < Xj-i and Sj < CN~ 1 t 2 . The 2" +1 
possible meshes are divided into subclasses of two types. On the meshes Mb 
with bi = the result follows from ([26]) . ([29]) and Lemma IT2l On the meshes 
Mb with bi = 1 the result follows from (|26|) . ([30]) and Lemma ITSl 
Thirdly, in the general case Xj S (T m ,T m +i) for 2 < m < n — 1, it follows 
that T m < Xj-i and Sj < C/V _1 T m +i. Then Mb is divided into subclasses of 
three types: Mj° = {M b : 6 X = • • • = 6 m = 0}, M£ = {M b : 6 r = 1, & r +i = 
• • • = 6 m = for some 1 < r < m - 1} and M^ = {M b : 6 m = 1}. On M£ 
the result follows from <j26j), (J29]) and Lemma [12] on M£ from (f26|) . ([29 ]) . p)]) 
and Lemma El on Af b " from ([H, ([30]) and Lemma US] 
Finally, for Xj G (r„, 1), r„ < a^-i and <5j < CN^ 1 . Then Mb is divided into 
subclasses of three types: Mg = {M b : &i = • • • = b n = 0}, M£ = {M b : 6 r = 
1, 6 r+ i = • •• = b n = for some 1 < r < n - 1} and M£ = {M b : b n = 1}. 
On M^ the result follows from ([26]), (J29]) and Lemma [H] on M£ from ((26j) , 
(J29]), (0Q]) and Lemma US] on M£ from ((30]) and Lemma US 

Now assume that Xj = t^. € Jb • Analogously to the proof of Theorem [1] 
the ratio R(u>i(Tk)) is introduced in order to facilitate the use of Lemma [TT] 
To complete the proof it suffices to establish in all cases that 



The required estimates of the denominator of R(iVi(Tk)) are ([43]) and (|44| . 
The numerator is bounded above using Lemmas Q2] and [TU The cases = 1 
and bk = are treated separately and the inequalities ([2"6"]l . ([2"T]) , ([2~5]) . ([3H]) 
and ([33]) are used systematically. 

Suppose first that bfc = 1, then there are four possible subcases: 



W - w|| < CN- 2 (\nN) 3 . 



(53) 



(54) 
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LemmaHl i < k, H k > h k , R{w r { Tk )) < C(f- + £M±±± ). 

H k < h k , R( Wl (r k )) < C(| + ( J9 3/2 ). 
Lemmadl * > k, H k > h k , i?K(r fc )) < C(l + v2EEH). ,; ' ;,! 

ff fe < flK(r fc )) < c(i + ^=). 

Secondly, if b k = 0, then b k -i = 1, because otherwise r k ^ and further- 
more H k < h k - There are two possible subcases: 

LemmaEl i < k - 1, H k < h k , < + 1). 

LemmalH » > fc - 1, -fffe < h k , R{wi{r k )) < C(l + §*). i ;>< " 

In all six subcases, because of the ordering of the it is clear that condition 
([Ml) is fulfilled. This concludes the proof. ■ 

The following theorem gives the required essentially second order parameter- 
uniform error estimate. 

Theorem 3. Let A(x) satisfy 0) and (0). Let u denote the exact solution 
from rip and U the discrete solution from \3J$ . Then 

||U - u|| < CN- 2 (\nN) 3 . (57) 

Proof. An application of the triangle inequality and the results of Theorems 
Q] and [2] leads immediately to the required result. ■ 
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